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Abstract 



We investigate the analyticity property of the partially resummed series ex- 
. pansion(PRSE) of the partition function for the eight-vertex model. Develop- 

ing a graphical technique, we have obtained a first few terms of the PRSE and 
found that these terms have a pole only at one point in the complex plane 
ON , of the coupling constant. This result supports the conjecture proposed by 



Guttmann and Enting concerning the "solvability" in statistical mechanical 
lattice models. 



Recently Guttmann and Enting introduced a new point of view into the study of the 



solvability of lattice models in the statistical mechanics, including models for combinatorial 
problems Their approach is based on a study of the connection between the solvability 
of lattice models by use of the inversion relation and the analyticity properties of their 
associated series expansion, which is called the partially resummed series expansion (PRSE). 
For exactly solvable models, the inversion relation for the partition function was derived 
using the star-triangle relation ||[| . Subsequently, it was shown to hold true for the partition 
function of (so far) unsolved models and other physical quantities as well [|]]. Thus the 
inversion relation can be used for clarifying the difference between exactly solved quantities 
and unsolved ones in statistical lattice models. 

The importance of the analyticity properties of the PRSE was first pointed out by 
Baxter ||. He considered the PRSE for the reduced free energy lnA(t 1; t 2 ) of the 



zero-field anisotropic Ising model on the square lattice, lnA(ti,t 2 ) = J2 RnitDtf 1 , and 
showed that the functions R n (t\) can be determined recursively from the inversion rela- 
tion, In A(ti, t 2 ) +ln A(£f , — 1 2 ) = ln(l — t|), an d the symmetry relation A(t 1; £ 2 ) = A(t 2 , ti), 
provided R n {t\) have singularities only at a single point in the t^-complex plane. This re- 
alization led to a series of investigations on the partial sums R n (t\) for unsolved quantities 
such as the partition function of the 2D Ising model in a field ||, non-critical 2D Potts 
model fU, zero-field 3D Ising model || or the susceptibility of the 2D Ising model Q. The 
results for these quantities indicate that infinitely many poles appear in partial sums R n in 
the limit n — > 00 in clear contrast with the solvable case. 

The new approach of Guttmann and Enting to the solvability of lattice models is as 
follows. They defined a given quantity to be solvable if its solution can be expressed with a 
D-finite function as a function of the coupling constant. D-finite functions are defined as 
functions which satisfy a linear ordinary differential equation of finite order with polynomial 
coefficients. The definition of the solvability was extended so that solvable quantities do not 
necessarily have a closed form solution. However, it is restrictive in the sense that those 
quantities which are expressible with the solution of a nonlinear differential equation are 
classified as "unsolvable" ones. 

Guttmann and Enting further put forward a criterion which determines whether a given 
quantity is solvable or not in the sense of D-finiteness. First, generalize the model such that 
it possesses two (or more) anisotropic coupling constants t\ and t 2 . Let the wanted quantity 
be /(ti,i 2 ). One can calculate exact coefficients of the series expansion of /, /(ii,t 2 ) = 
^a mn t™t2 (to a considerably high order with sophisticated techniques). Then, performing 

mn 

the infinite sum with respect to, say, t\ first, one obtains the PRSE, f{t 1 ,t 2 ) = ^ Rn{ti)t\. 

n 

Let us assume that the poles in R n (ti) all lie on the unit circle in the ti-complex plane (this 
is not essential, however). Their conjecture on the criterion is stated as follows: (i) If the 
number of poles of R n remains finite as n —>■ 00, then /(ti,t 2 ) is a D-finite function of ti 
with t 2 being fixed, (ii) If there appear infinitely many poles densely distributed along the 
unit circle as n — > 00, then /(ti,t 2 ) cannot be expressed in terms of D-finite functions only. 
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It is illuminating to look at the spontaneous magnetization M of the 2D Ising model as an 
example of the case (i). The closed form solution is given as M(u, v) = (1 — 16uv /(l — u) 2 (l — 
w) 2 ) 1 / 8 in the low temperature phase, where u — (1 — ti)/(l + ti) and v — (1 — t 2 )/(l + h)- 
One can calculate a first few terms of the Taylor expansion M(u,v) = J2 n >o Rn{u)v n by 
directly taking derivatives of the above expression. But this would soon become undoable. 
A concise way to calculate the coefficients is to make use of the differential equation satisfied 
by M(u,v). By extending the argument in |§ one arrives at 

d 

(1 - u) [(1 - vfu 2 - 2(1 + 6v + v 2 )u + (1 - vf] jr-Miu, v) + 2v{l + u)M(u, v) = 0. (1) 

CJ Lb 

This shows that the spontaneous magnetization is D-finite. It is then determined implicitly 
from the recursion relation for R n 

(n + 3)(1 - u) 2 R n+3 (u) - [3(n + 2) + 2(5n + 9)u + 3(n + 2)u 2 }R n+2 (u) 
+[3(n + 1) + 2(5n + 6)u + 3(n + l)u 2 ]R n+1 (u) - n(l - u) 2 R n (u) = 0. (2) 

From this recursion relation, it is clear that the partial sums R n have poles only at u = 1. 
We should note that the recursion we have derived above from D-finiteness is different from 
that discussed by Baxter H and mentioned earlier. In the latter case, to calculate the term 
R n one needs all Rk prior to R n . By contrast, in the former case, only a fixed finite number 
of terms is sufficient to obtain the term R n no matter how large n is. 

Examples of the case (ii) are the zero-field susceptibility of the 2D anisotropic Ising model 
and generating functions for self-avoiding polygons on the square and hexagonal lattices. 
Guttmann and Enting argued using the finite-lattice method || that infinitely many poles 
will appear in partial sums for those quantities in the limit n — > oo [|TJ . 

The few analyses mentioned above already indicate that Guttmann and Enting's con- 
jecture unravels a new aspect of exactly solvable models, the connection of the solvability 
to the distribution of poles and D-finiteness of the solution. The conjecture deserves a fur- 
ther and detailed study. In particular, it is very important to examine whether Guttmann 
and Enting's conjecture has general applicability to a broad class of statistical mechanical 
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systems. The eight-vertex model is particularly suited to this end. It can be regarded as 
a nontrivial generalization of the 2D Ising model. It has three coupling constants and can 
be compared with the previous analyses of the zero-field 3D Ising model and the 2D Ising 
model in a field. 

We adopt the spin-type model description of the eight-vertex model. The spin variable 
a placed on each site of the square lattice consisting of iV faces(Fig.l) takes the values ±1. 
The local Boltzmann weight is defined on each face / in the lattice and depend on the spin 
configuration around the face {a} / = (a, b, c, d) as 

w({a} f ) = Rexp(Kac + Lbd + Mabcd). (3) 

The partition function is given by 

Z N (K,L,M) = X)n w ({°'}/)> ( 4 ) 

where the summation is over all spin configurations on the lattice sites and the product is 
over all faces in the lattice. The lattice is assumed to have the periodic boundaries. This 
model can be regarded as a pair of Ising models with anisotropic coupling constants K and 
L, each living on one of the sublattices and coupled by the four-spin coupling M. 
The local Boltzmann weight can be rearranged as 

w({a} f) = p(l + ziac + z 2 bd + z 3 abcd), (5) 

where 

tl + #2^3 ^2 + ^3^1 ^3 + tlt 2 

1 + M 2 t 3 I + M2V I + M2V 

p = i?(cosh K cosh L cosh M + sinh K sinh L sinh M) , (7) 

with ti = tanhfT, £2 = tanhL, £3 = tanhM. z\, z%, £3 are high-temperature variables and p 
is the normalization factor. Let us represent each term in RHS of ([5]) graphically by bonds 
as shown in Fig.2. We call the bonds shown in (b), (c), (d) of Fig.2 z±-, Z2-, ^-bonds. 

Substituting (^) into (|) and expanding the products yields all possible bond configura- 
tions over N faces in the lattice. The summation over spin configurations leaves only those 
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in which an even number of bonds are incident to each lattice site, i.e., closed graphs. The 
reduced partition function, A N (zi, z 2 , z 3 ) = (2p)~ N Z N (K, L, M) , also plays the role of the 
generating function for the enumeration of closed graphs in the lattice. They are classified 
in terms of the numbers of the Zi~, Z2-, Z3-bonds, l,m,n. Then can be written as 

00 00 00 

A N ( Zl , Z 2 , ^ 3 ) = E E E Plmn{N)zl l zl m zl\ (8) 
1=1 m=l n=0 

where Pi mn {N) is the number of the closed graphs and is a polynomial in N. Define 
the reduced partition function per face in the thermodynamic limit by A{z\,z 2 ,z 3 ) = 
lim (Ajv(zi, z 2 , z 3 )) l ^ N . Then the expansion of A(zi, z 2 , z 3 ) in the connected closed graphs 
is given by 

00 00 00 

In A( Zl , z 2 , z 3 ) = E E E a lmn zfz 2 2 m z 2 3 n , (9) 

1=1 m=l n=0 

where a\ mn is the coefficient of N in the polynomial Pi mn fllQ| . 

Our aim is to see whether or not lnA(z 1 ,z 2 ,z 3 ) possesses the property conjectured by 
Guttmann and Enting. We begin by looking at the inversion relation for the eight-vertex 
model 0, which reads in the variables zi, z 2 , z% as 

In A(zi , z 2 ,z 3 ) + In A M ^-^ J Zl , -z 2 , - ^ 1 ~ ^ =ln(l-2; 2 2 ). (10) 

Note that if one sets the four-spin coupling M to zero, i.e. z 3 = Z\Z 2 , (|I0|) reduces to the 
inversion relation for the anisotropic Ising model. Since (|T0|) involves the inversion of two 
variables, Z\ and z 3 , one has to take the double infinite sum with respect to both to compare 
the series (^) with (|T0|). In practice this turns out to be a hard task. Fortunately, one can 
circumvent the difficulty by a slight change of variables as V\ — z\,v 2 — z 2 and v 3 = z 3 jz\. 
Then ( |T0D becomes 

lnA{v u v 2 ,v 3 ) +lnA(^^0v^\-v 2 ,-v 3 S] j =ln(l-w 2 2 ). (11) 

In this case, only v% is inverted and one has only to take the partial resummation with 
respect to vi, which we write as 
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hlA{v 1 ,V 2 ,V 3 ) = J2 Y. R mn{vl)vl m vl n . (12) 

771=1 71=0 

The corresponding expansion of the second term in LHS of ([11]) can be obtained by 
computing the Taylor expansion of R mn ({l — ff) 2 /[vi(l — v 2 )] 2 ) with respect to v\ and v\. 
Collecting the terms containing the same powers in v\ and v 2 , one has the expansion, 

k a((J— 4V\ -V2, -v 3 ) = E E g m n(^)^ 3 2n - (is) 

\\l — v 3 ; j m=1 n=Q 

The inversion relation (|ll|) is now written as a perturbative relation, 

Rmn{vl) + Qmniyl) = — <W (14) 

m 

Note that Q mn depends only on R^i with k + I no greater than m + n. The relation (|14"D 
gives constraints on the partial sums R m n(v 2 ) and is used in the later analysis of the PRSE. 

One might think that it is immediate to obtain the PRSE (|i~2|) for the eight- vertex model, 
since we have the exact solution at hand. However, the expression of the exact solution in 
terms of the elliptic parameters is very implicit. So far we have not been able to derive the 
PRSE from the exact solution and to find out the location of poles in the partial sums in 
v 2 . 

Instead we employ a graphical approach to calculate R mn {v 2 ) in the series expansion 
(|12|). The closed graphs with 2m Z2-bonds, 2n z 3 -bonds and an arbitrary number of Z\- 
bonds contribute to R mn . Because of the presence of the four-spin coupling constant, there 
appear vastly many types of closed graphs even at low orders. We have developed a scheme 
to classify the graphs by introducing the notion of minimal graph, the graph with the least 
number of 24-bonds among those in one class. The derivation of our method is lenghty and 
will be presented in a separate paper. Here we report the results of our analysis of a first 
few terms of the expansion: 

m = l,n = The minimal graph is the square(Fig.3(a)). The remaining graphs are 
obtained by successively adding a pair of Zi-bonds vertically, which we call u^-extension. 
The weights of the i^-extended graphs in (a) of Fig. 3 are v\v\, vfv 2 , v\v\ and so on. Their 
contribution to R w sums up to v\/(l — v\). 
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m = l,n = 1 There are two classes of graphs with two Z2-bonds and two 23-bonds 
which contribute to Ru. One of them is shown in (b) of Fig. 3. Its minimal graph can 
be extended in three ways. It has two t> 2 -extensions, extending either of the two squares 
independently. The third is a simultaneous extension of the two squares (see the third 
graph of Fig.3 (b)), which we call v ^-extension. The weights of the graphs in (b) of Fig. 3 
are thus vfv%v%, vfv%v%, vfv^v 2 and so on. The contribution from these graphs sums up to 
2vf/(l — vf) 2 (l — vf), where the factor 2 takes account of embeddings of graphs and their 
reflection. The contribution from the other class to Rn is given by 2vf/(l — vj) 2 (l — vf). 

m = 2, n = There are six classes of graph with four Z2-bonds and zero 23-bonds which 
contribute to i?2o- Five of the corresponding minimal graphs consist of connected graphs and 
the other of a disconnected graph with two squares. These graphs have only v\ -extensions. 

m = 1, n = 2 There are twelve classes of graphs which contribute to R\2- In each class 
the minimal graph can be extended in many different ways. Some of the twelve minimal 
graphs are shown in Fig.3(c)-(e). The graphs of the class (c) are obtained from the minimal 
graph (c) by two ways of v ^-extension and three ways of t^-extension. Their contributions 
yield 2vf/(l — v 2 ) 2 (l — vf) 3 . The contributions of the graphs of the class (d) and (e) can be 
found by similar arguments: 4vf/(l — v 2 ) 2 (1 — vf) 3 and 2v\°/(l — v 2 ) 2 (l — vf) 3 respectively. 
The sum of the contributions from the twelve classes turns out to take a surprisingly simple 
form, 3vf(l + vf)/(l — vf) 5 , which is free from the pole at v\ = —1. 

To summarize the above analysis, we have obtained the following expressions for the first 
four terms of the PRSE: 



l _ v 2> (l-„2)3' 

_ v 2 (2-5v 2 + vf) _ 3v*(l + vf) 

R ™y v i) - — (1-^2)3 — ' - (i- v f)5 ■ ^ 

The corresponding terms of Q mn (vf) can be easily obtained from these results as explained 
previously. 

There are a few important points to be noted. First, we have checked that R mn {vf) 
and Q m n( v i) satisfy the perturbative form of the inversion relation (0). This is a highly 



nontrivial check that the result fll5|) is correct. The second point is concerned with the pole 
structure of R mn (vf). As remarked above, Rio(vf) and R2o(vf) have a pole only at v\ — 1 
for a graphical reason. As for Rn(vf), both of the contributions from the two classes have 
poles at vf — 1 and v\ = — 1. Curiously, the extra pole at v f = —1 cancels out in the sum 
of two contributions, leaving the pole term 1/(1 — vf) 3 . It is even more miraculous that, as 
noted above, the same kind of cancellation of the poles at v\ = —1 takes place for Ri2(vf), 
leaving the pole term 1/(1 — vf) 5 . 

We think that the absence of extra poles at v \ = — 1 in the first four terms of the PRSE 
is not a pure accident but that it has a raison d'etre. It should be counted as a strong 
indication that poles in R mn (vf) will not appear at other values than v\ = 1 for higher m 
and n, supporting the conjecture of Guttmann and Enting. 

It is illuminating to compare the present result on the PRSE for the eight-vertex model 
with that for the 3D Ising model ||. The latter model also has three coupling constants. The 
inversion relation and the PRSE in this model read as In A(t 1; t 2 , t 3 ) + lnA(£f , —t 2 , —t 3 ) = 
ln(l — t|) and In A(t 1 , t 2 , t 3 ) = R mn {t^)tl m t^ n . Nontrivial results start to appear from 

mn 

Ruitf) (Rm,o(vf) = R 0m (vf) coincide with R m {t\) of the 2D Ising model, Ru{t\) has no 
poles other than at t\ — 1 for a straightforward graphical reason). Indeed, Ruitl) was 
found to have poles of the form 1/ (1 — tf) 5 (1 + tf) . An extra pole at t\ = — 1 already appears 
at the third order (m + n = 3) term of R mn (t\). It is a sign that an infinite number of poles 
will appear in the limit m, n — * 00. We note a clear difference in the pole structure of R mn (t) 
between the eight-vertex model and the 3D Ising model even at low orders. Probably this 
difference reflects the fact that the former model is solvable while the latter is not. 

Finally, we postulate the form of the partial sum R mn (vf) at all orders by extending our 
results at low orders and examine the role of the inversion relation regarding the exactly 
solvable partition function in the approach of the PRSE. Let us assume that the partial 
sums of the eight- vertex reduced partition function have only (multiple) poles at vf — 1. By 
a graphical consideration, one can guess the general form of the partial sums: R mn {vf) = 



S 



2(m+n)-l 

P mn (vf)/(1 - v f) 2 ( m + n )- 1 j where P m n{vf) is a polynomial in vf , P mn (vf) = c mnkvf 

fc=0 

Then P mn (vf) can be determined recursively from the inversion and symmetry relations 
as follows. Suppose that Rm with k + I < L — 1 are known. Then from the symmetry, 
In A(t> l5 1> 2 , ^3) = In A(t> 2 , v±, Viv 3 /v 2 ), we have the system of linear equations, 

L-i 1 2(L - 1) + r - k \ 2r-i / 2(r - 1) + L - k \ 

fc=o V 2(L - 1) / k=o V 2(r - 1) / 



where < r < L — 1 and c mn fc = if m < 1. Since the RHS of ( |16|) is written with known 
coefficients only, we have a half of the coefficients {c mn k :m + n = L,0<k<L — 1}. From 
the inversion relation (0) we have 



Rmn{vl) = -R mn (v 1 2 ) + known terms, (17) 

which determine the other half of the coefficients {c mn k : m + n = L,L < k <2L — 1}. Thus 
the solution for the eight-vertex partition function is reproduced by extending the method 
which was first used for the 2D Ising model by Baxter [|[ . The first few terms of this solution 
agree with (0). Moreover, it reduces to the PRSE for the Onsager solution in the Ising 
limit. This remarkable coincidence further supports our argument on the pole structure of 
R mn for the eight- vertex partition function. 

The authors thank Makoto Katori for bringing the paper of Guttmann and Enting to 
their attention. They are indebted to Atsuo Kuniba for discussion at various stages of the 
present work, and to Satoru Odake for valuable comments on the manuscript. They are also 
grateful to Tony Guttmann for elucidating conversations, in particular on the differential 
equations to be obeyed by solvable quantities. 
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FIGURES 



FIG. 1. The eight-vertex model. 



FIG. 2. Bond configurations on a face. 



FIG. 3. Minimal graphs and their extensions. 
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